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ABSTRACT 

This paper provides the theoretical framework for the development of a 
pseudo-ciosed loop solar tracking system incorporating a microcontroller 
for a fiat payload, where the outputs of the system (the slope and the 
azimuth angles of the payload) are observed and the set-points for the 
system are to be optimally computed by the microcontroller, it also 
presents, from literature, the equations and models solar angles and 
incident angles for the development of the optimum set point path for the 
tracking, and gives the outline of the signal flow plan, as well as the 
corresponding flowchart, for the tracking system's components interaction. 
The concept would provide a cheaper and simpler alternative to dual axis 
solar tracking incorporating solar sensors. The models presented here are 
adaptable to any geographical location in Nigeria and globally. 

© 2019 Faculty of Engineering, University of Maiduguri, Nigeria. All rights reserved. 

1.0 Introduction 

Conventional energy sources like oil and coal still constitute more than 80% of the global 
energy mix (IEA, 2018). Yet apart from being finite, the conventional energy sources pollute and 
warm the environment. According to a report by (NASA, 2018), about 97% of actively publishing 
climate scientists agree that the earth's climate is warming due to human activities, especially 
the use of fossil fuels. These reasons have given impetus to the exploration and use of non- 
conventional, clean, and renewable energy sources such as solar energy and wind energy. Solar 
energy is the direct utilization of the sun's radiant energy. 

A major component of all solar energy applications is the payload - the part that absorbs the 
solar energy directly and converts it to a useful form, or transmits it to another media. In helio¬ 
electrical applications, the payload is the solar photo-voltaic panel, and it is a solar collector in 
helio-thermal applications. Solar collectors absorb solar energy and convert it into heat which is 
then transferred to a stream of liquid or gas. To get the most from solar payloads, they must be 
pointed to the directions that capture the most solar energy. In general, results show that solar 
tracking can increase the yearly energy capture by up to 40% compared to the case in which the 
payload is fixed (Bari, 2000, Lee et al„ 2009, Mohd et al., 2001). 

Solar tracking can be manual or automatic. In automatic solar tracking, the payload tracks the 
sun at all times with very little human intervention (that is the angle of incidence of the sun with 
respect to the receiving surface of the payload equal to zero at specified times). In the strict 
sense, there are two categories: open loop tracking and closed loop tracking. Between these 
two categories, there can be another - a pseudo closed loop solar tracking. All these are 
explained as follows. 
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1.1 Open loop solar tracking 


In open-loop tracking, data of the desired payload's altitude and surface azimuth angles at 
specified times are computed based on some sun-earth geometry or models, or solar irradiation 
algorithms. The data is fed into, or computed by, an appropriately coded microcontroller which, 
using information from a real time clock sends signal to an actuating mechanism to manipulate 
the altitude and surface azimuth angles of the payload appropriately. Figure 1 gives a schematic 
of a general open-loop solar tracking. It can be seen from the figure that an open-loop solar 
tracking system only manipulates the altitude and/or surface azimuth angles of the payload, it 
does not observe the actual payload's incidence angle using any form of sensor or feedback 
mechanism, so it cannot correct for any errors in the actual incidence angle. For this reason, it is 
poor in disturbances rejection. 



Actual altitude/ 
azimuth angle 


Figure 1: Schematic of open loop solar tracking 


Models and algorithms for computing the desired paths for open-loop solar tracking systems 
exist in literature. Some of these modes and algorithms require basic input parameters such as 
date, time (local standard time (LST), time zone (TZ)), longitude, latitude and surface azimuth 
rotation of a location. The very complex and precise algorithms are contained in Blanco-Muriel 
et al. (2001), Michalsky (1988), Grena (2008), Meeus (1998) and Reda and Andreas (2004). Actual 
development or design and construction of open loop solar trackers using a form of the 
formulas and algorithms already mentioned have also been reported. Among them is Al-Naima 
and Yaghobian (1990), who designed and constructed a dual axis sun tracker in which the 
controller is a microprocessor. Abdallah and Nijmeh (2004) used PLC controllers in the dual sun 
tracker they developed. 

1.2 Closed loop solar tracking 

Closed-loop tracking systems are based on feedback control principles. In these systems, the 
payload's incident angle is sensed indirectly through a system of light detecting or measuring 
devices such as light detecting resistors (LDR), pyrheliometer and four quadrant photo detectors. 
Using feedback mechanisms, the error signal from the sensed incident angles are sent to a 
controller, which carries the necessary control action and send to the actuator. A form of closed- 
loop tracking system (with provision for alternate set point input) is shown in Figure 2. Because 
of the light sensor and controllers required in their operations, closed loop solar tracking 
systems are generally more complicated and expensive than the open-loop solar tracking 
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Different types of closed loop solar tracking systems have been developed. Kalogirou (1996), 
Barsoum (2011) and Aiuchi et al. (2004) employed light-dependent resistors (LDRs) in their one- 
axis closed loop solar tracking systems. Jiang and Cao (2003) used an emulated sunflower based 
on a spherical four-quadrant photoelectric sensor for sensing the incident angle of their solar 
tracking systems. Roth et al. (2004) and Georgiev et al. (2004) designed and constructed a sun 
tracking system in which a pyrheliometer was used to measure the direct solar radiation, a four- 
quadrant photo detector to sense the sun's position. Xinhonz et al. (2007) and Alata et al. (2005) 
used fuzzy logic algorithms in their controllers. 

The objective of this paper is to develop a tracking system that is a bridge between the open 
loop system and closed loop system described above. Such developed tracking system is 
termed pseudo-closed loop solar tracking system and the process is presented in the following 
Section. 

2.0 Materials and Methods 

The main material for this research is a computing facility with MATLAB® software. The 
methods include conceptualising the pseudo-closed loop solar tracking system, presenting and 
simulating the solar energy equations required for developing the optimum set points for the 
pseudo tracking system, and developing the procedure required for implementing the pseudo 
closed loop solar tracking system. 

2.1. The pseudo-closed loop solar tracker concept 

The proposed pseudo-closed loop solar tracker will have all the features of an open loop system 
but would also incorporate a feedback system to observe the altitude and surface azimuth 
angles of the payload. It is a pseudo-closed loop tracker in the sense that the feedback would 
not observe the actual incidence angle of the sun. The development of the concept of the 
tracker, and its components, as well as the development of the equations to be used for 
obtaining its optimum set-points, and its signal flow plan, are presented in the following Sub- 
Sections. 

Figure 3 is a schematic of the pseudo-closed loop solar tracker. Compared to a closed loop 
system of the type shown in Figure 2, the pseudo-closed loop system would use a smaller 
number of sensors/transducers, and would avoid the use of controllers like PID, Fuzzy logic and 
the like, that may introduce additional complications in design and implementation. It has as its 
major components a microcontroller, a real time clock, angular position sensors and actuators. 
The microcontroller uses signal from the real time clock to determine which value from the 
optimum altitude and azimuth data set to be used as error signal for the actuator at specified 
time interval. The optimum dataset may be supplied externally, as indicated in Figure 3, or the 
codes for generating it may be stored in the microcontroller and the optimum dataset 
computed 



Figure 3: Schematic of pseudo-closed loop solar tracking system 
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This concept does not correct for any errors in the actual incidence angles: this will depend on 
the accuracy of the algorithm used to obtain the optimum set point paths, and the system set¬ 
up and calibration during installation. In addition, having optimized set-point calculation 
incorporated in a pseudo-closed loop system is an avenue to experiment on the correctness of 
existing solar angle models and algorithms. 

2.2 Developing the Optimum set point path for the pseudo-closed loop solar tracker 
The main aim of solar tracking is always to try and keep the angle of incidence of solar radiation 
on a payload at or close to zero . Optimum set point path is about being able to specify the 
degree of vertical and/or horizontal rotation of the payload necessary to make the angle equal 
to zero at specified times. For this a form of solar geometry formulas or algorithms for 
estimating the angle of incidence and the position of the sun, is required. 

For many solar energy applications, tracking accuracy of between 0.1° and 0.25° is sufficient 
given that the diameter of the solar disk subtends an angle of 0.52° (at aphelion, when the sun 
is farthest) to 0.54° (at perihelion, when the sun is closest) in the sky (Myer, 2013). For these 
situations, solar tracking equations or algorithms used need not be overly complex. The basic 
equations for estimating the angle of incidence as contained in Duffie and Beckman (1991), 
Tiwari (2002), Myer (2013), and some solar energy books as described in this section will suffice 
for those situations. 

2.2.1 Angle of incidence (00 

For a flat payload, the angle of incidence (9;) is the angle between the beam radiation and the 
normal to the surface of the payload (Duffie and Bechman, 1991). This angle is a function of 
three sun-earth angles - the latitude angle (4>), the declination angle (5), the hour angle (u>), 
and two other angles which define the placement of the payload on a surface - its slope angle 
(P) and its surface azimuth angle (y). The geometrical relationships between these angles are 
shown in Figure 4. The angle of incidence is defined by its cosine (Duffie and Beckman, 1991) as: 

cos0i = sinSsincpcosP + cosScoscpcosPcosoi + cosSsincpsinpcosycosco m 

4 - cosSsirLpsinysinco — sinScoscpsinpcosy ‘ 1 ' 

The latitude angle, together with the longitude angle, uniquely defines a location on the earth's 
surface. The declination angle and the hour angle are defined as follows: 

a) Declination angle (6) 

The declination angle is the angle between the equatorial plane and a line joining the centres of 
the Sun and the Earth, and this angle gives rise to variation in the length of the day (sunshine 
hours) throughout the year. The equation for calculating the declination angle for a non-leap 
year is given (Duffie and Beckman, 1991) as: 

6=23 M(^) (284+n) ] (2) 

where n is the number that represents the day of the year, for example, n = 32 on the 1st of 
February. 
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Figure 4: Relationship between Azimuth, Zenith, Altitude and Incidence angles 
b) The Hour angle (w) 

The hour angle is the angle of the sun from Solar Noon. For solar radiation calculations, the 
solar time is used instead of the local standard time. When the solar time (ST) for the location in 
hours is known, the hour angle can be calculated from (Tiwari, 2002): 

co = 15 (ST-12) (3) 

While the local standard time (LST) is always known, the solar time (ST) for a location is 
obtained from calculation as (Iqbal, 1983): 

ST = LST - 4(i)j st - i|/ a ) + ET (4) 

where, 

i|/ st = the standard longitude (time zone), in degrees, for the location. 
i|/ a = the actual longitude angle (in degrees) for the location. 

LST = the local standard time (in hours) 

and ET is the equation of time in hours given as (Spencer, 1971): 

ET = plr) (0.000075 + 0.001868cosB-0.032077sinB-0.014615cos2B- (5) 

0.04089sin2B) ' 

where B is defined as: 

B = 360(n — l)/365 (6) 

2.2.2 The position of the sun 

The position of the sun with respect to a location on the earth's surface can be defined using 
the altitude angle (a), the Zenith angle (0 Z ) as well as the solar azimuth angle (y s ). The 
relationships between these angles are defined as follows: 

a) Altitude angle and Zenith angle 

The zenith angle is the angle between the sun's ray and the perpendicular line to a horizontal 
plane at the location (Figure 7). The altitude angle (a) is the complement of the zenith angle. 
The altitude angle (and the zenith angle) can be calculated from the expression (Myer, 2013): 

sina = cos0 z = sinSsincf) + cosScosaicoscf) (7) 

b) The solar azimuth angle (y s ) 

In the northern hemisphere, the solar azimuth angle is the angle between the projection on a 
horizontal plane, of a line between the sun and a point (say P), and the line due South (Figure 4). 
In the southern hemisphere, it is defined as the angle between the projection and the line due 
north. By convention, the azimuth angle is taken as negative East of South and positive West of 
South for northern hemisphere (it is taken as negative East of North and positive West of North 
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for southern hemisphere). This means that y s is negative in the morning, and positive in the 
afternoon, and therefore, the sign of y s should match that of the hour angle o>. The solar 
azimuth angle can be calculated from the expression (Myer, 2013): 

cosy s = (sinasincf) - sinS)/(cos5cos4>) (8a) 

or 

cosYs = (cosSsincfjcosa) - sin5cos4>)/cosa = (cosSsincpcosa) - sin5cos4>)/sin0 z (8b) 
Taking the inverse cosine to obtain y s in Equation (8) will not always come with the correct sign. 
So, it must be ensured that y s takes the sign of the hour angle oo as explained before. If 5 > 0, 
the sun can be north of the east-west line. 

In addition to the angles mentioned in the above equations, it is important to know when it is 
sunrise or sunset, and calculate its accompanying angles, so that the actuator may be instructed 
about periods of idleness or activity. For this the following angles and hours are defined. 


c) Hour angle at sunrise (o) sr ) and sunset (o> ss ) 

At sunrise, the zenith angle is 90o (altitude angle equals to Oo) on a horizontal plane, and this 
information is used in Equation 7 to obtain the expression for the hour angles at sunrise and 
sunset. Thus Sunrise (or sunset) hour angle can be calculated from (Tiwari, 2002): 

oo s = cos _1 (— tancptanS) (9) 


The sunrise (or sunset) angle can also be evaluated using the following equations (Myer, 2013): 
cos« z - sindsincp 


cosScoscf) 


( 10 ) 


cosasiny 

cos5 


( 11 ) 


The hour angles at sunrise ( oo sr ) and sunset ( oo ss ) are the same numerically (i.e. 
|oL» sr | = |co sr | = oo s ), however the sunrise hour angle is negative and the sunset hour angle is 
positive. 


d) Number of daylight hours (N) 

The total angle between sunrise and sunset, in degrees is thus equal to 2o> s . This angle between 
sunrise and sunset (2o) s ) can be used to find the number of daylight hours (N) for a particular 
day by using the expression (Tiwari, 2002): 


There are always 4380 hours of daylight per year (non-leap year) everywhere on the globe. 
When (- tan5tanc|)) > 1, there is no sunrise (i.e. 24 hours of darkness). This implies that oo sr = 0. 
Also the solar times at sunrise and sunset for this condition are both equal to 12. When 
(-tanStancf)) <— 1, there is no sunset (i.e. 24 hours of daylight). This implies that a> sr =-180 
and co ss = 180. The solar times at sunrise and sunset for this condition are respectively equal to 
0 and 24. 


2.2.3 Optimum set point path for solar Tracking. 

Equation (1) shows that in order to make a flat payload track the sun, only its surface orientation 
(P) and surface azimuth angle (y) are open to adjustments. The others (the latitude, the 
declination and the hour angles) are automatically determined from knowledge of the location 
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and the date and time. Therefore, the cosine of the angle of incidence may be expressed as a 
function of two variables by re-writing Equation (1) as: 

cosBj = Cicosp + C2sinpcosy + C3sin(Bsiny (13) 

where 

Ci = sinSsincf) + cosScos 4 >coso) 
c 2 = cos5sin4>cos(jL> - sinScoscf) 

C 3 = cosSsinco 

From Equation (13), the pair of the values of (B and y which makes the cosine of the angle of 
incidence equal to unity, or the angle of incidence equal to zero, at any location and for 
specified date and time, is that which makes the payload perfectly tracks the sun at that time. So, 
sun tracking is about finding the optimum values of this pair at desired times. 

Equation (13) can further be reduced to one of a function of a single variable if, at any location 
and for specified date and time, the payload is oriented such that its surface azimuth angle (y) is 
numerically equal to the solar azimuth angle (y s ), as can be seen from Figures 4 and 5. So 
Equation (13) may be re-written for case the of perfectly tracked sun as: 

c 6 cosp + Cysinp - 1 = 0 (14) 

where 

C 6 = sinSsincf) + cosScos4>cosa> 

c 7 = cos5sin4>cos(jL>cosy s + cos5sina)siny s - sin5cos4>cosy s 

With reference to Figure 3, a microcontroller in a pseudo closed loop control can be coded to 
calculate the desired inclination angle ((B) and azimuth angle (y) a location at regular time 
intervals in a day, or even as once in a day, once in a month or once a year calculation. Then this 
optimum data set would be sent to the actuator for the desired altitude and azimuth alignment. 
The feedback path in the figure (with its sensors or transducers) is optional, just to provide 
either historical data and/or to provide information about the correctness of the control system. 
Optimum value of [B for known y s , at any location and period, can be evaluated from Equation 
(14) by closed-domain method for solving nonlinear equations such as bisection method or 
false position method. The MATLAB® software has a function fminbnd for solving single 
variable closed-domain problems. Using MATLAB® software, the codes for evaluating Equation 
(13), and all the associated equations from Equations 1 to 12, are implemented and the trends 
of the angle of incidence for a range of surface orientation ((B) and surface azimuth angle (y), for 
Kano, Nigeria (11.9683 °N, 8.4261 °E) at 14:27 hour on 19-OCT-2018, as an example, computed. 

2.2.4 Development of the signal flow plan for the pseudo-closed loop solar tracking system 
A signal flow plan, which describes an option of calculating the optimum angles for sun tracking 
as once a day, was developed. A flowchart for implementing such code for a microcontroller 
was also developed. The option of once a day calculation was preferred over periodic 
calculations within a day to reduce the computational load on the microcontroller. It was also 
preferred to the other options (once a month or year calculations) for data storage issues. The 
proposed pseudo-closed loop solar tracking system would operate by having its set-points 
(optimum azimuth and slope angle) being automatically adjusted. The development of the 
signal flow plan for the system is such that it should be fully operational for all days of the year, 
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if it is powered on. The complete signal flow chart is given in Figure 5, and the steps in the 
signal flow plan are outlined as follows: 



Calculate sunrise 
and sunset times 



Calculate 
Optimum 8 and y 


T 


Send 

Optimum data to 
internal storage 



Make current 





Figure 5: Flow chart for the pseudo-closed loop sun tracking system 
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At power on, the system reads the latitude and longitude data for the location. It then makes 
the set-points for the system (azimuth and slope angles) equal to zero. 

The system then reads the current date from the attached real time clock 

The system then calculates the sunrise and sunset times for the day, using Equations 2,3 and 9. 

The system then reads the current time 

The system then checks if the current time is within sunrise and sunset times calculated in step 3. 
If the current time IS NOT within sunrise and sunset times, the system calculates the necessary 
delay time td os based on the current time and the sunrise and sunset values calculated. The 
system is then delayed by td os , then returns to step 3 to attempt another date reading. 

If the current time IS within sunrise and sunset times, the system calculates the optimum 
azimuth and slope angles for specified time intervals (say every 4 minutes) between sunrise and 
sunset. It uses Equations 2-6, 8, 14 

The system sends the optimum set-point data to the internal storage of the microcontroller. 

The system reads the current time again (to compensate for any time lapse since the first 
reading) 

The system checks if the end of day has been reached (i.e. if current time is beyond the sunset 
time value already calculated. 

If IT IS end of day, the system calculates the necessary delay time td ed based on the current time 
and the sunset values calculated. The system is then delayed by td ed , then returns to step 3 to 
attempt another date reading. 

If IT IS NOT end of day, the system searches for a match between the current time and the 
stored times in the optimum data table of the microcontroller. 

If equality IS NOT found in step 11, the system delays by 30 seconds and return to step 8 to 
read the current time. 

If equality IS found, the system sets the corresponding slope and azimuth data that is a match 
for the data (in step 11) as the set point for the system. 

The system delays by 2 minutes (assuming the optimum times were calculated at 4-minute 
intervals) before returning to step 9. 

3. Results and Discussions 

As an example, the plots of the trends of the incident angles (0;) for various azimuth angles (y) 
and inclination angles ((3) are plotted as shown in Figure 6 for Kano, Nigeria (11.9683 °N, 8.4261 
°E) at 14:27 hour on 19-OCT-2018. While Figure 6(a) shows the 3-D plot, Figure 6b shows its 
contour plot. The contour plot (with the colour bar at the right) shows that at this time, the 
optimum pair of (3 and y lie in the darkest part of the plot (the dark blue region is the region of 
lowest angle of incidence), which evaluates to between 22° to 55° for (3 and 30° to 75° for y. 

While Figure 6 shows the different combinations of azimuth angles (y) and inclination angles ((3) 
to give possible incident angles (0;) at an instant (a given time), Figure 7 shows the plots of 
optimum pairs of (3 and y for a day. For the date specified for Kano, Nigeria on 19-OCT-2018, 
the optimum pair of (3 and y that will make a flat payload track the sun from sunrise to sunset 
are given in the figure. The solar time is shown on the horizontal axis while the azimuth and 
inclination angles are shown on the vertical axes. For example, the figure shows that the 
optimum pair of [3 and y at 10.00 solar time on this date is about 55° for [3 and -40° for y (the 
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points where p and y curves intersect the 10.00 solar time line). Likewise on this date, the 
optimum pair of p and y at 12.00 solar time is about Oo for p and 35° for y, and at 14.00 solar 
time, the optimum pair is about 45° for p and 40° for y. 

4.0 Conclusions 

A concept of pseudo-closed loop solar tracking system for a flat payload, where the set-points 
for the system are optimally computed by a microcontroller, and the outputs of the system (the 
slope and the azimuth angles of the payload) observed, has been proposed. The models for the 
development of the optimum set point path for the tracking system for use within a 
microcontroller are also presented. Example of using the models to calculate optimum payload 
orientation angle and azimuth angle for a given day was presented. The outline of the signal 
flow plan for the tracking system's component interaction, and the corresponding flowchart, 
were also developed and presented. The pseudo-closed loop solar tracking system offers a 
viable and cheaper alternative to a full closed loop solar tracking system. 



Figure 6: 3-D plot (a) and contour plot (b) of angles of incidence for ranges of inclination angle 
(p) and azimuth angle (y) for Kano, Nigeria (11.9683 °N, 8.4261 °E) at 14:27 hour on 19-OCT- 
2018 



Figure 7: Plots of optimum pairs of inclination angle (P) and azimuth angle (y) from sunrise to 
sunset for Kano, Nigeria (11.9683 °N, 8.4261 °E) on 19-OCT-2018 
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